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Mult i— frequency spectral analysis of extragalactic radio 
sources in the 33— GHz VSA catalogue: sources with 
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ABSTRACT 

We present an analysis of the spectral properties of the extragalactic radio sources 
in the nearly-complete VSA sample at 33 GHz. Data from different surveys are used 
to study source spectra between 1.4 and 33 GHz. We find that, in general, spectra can 
not be well described by a single power law in the range of frequencies considered. In 
particular, most of the VSA sources that are steep between 1.4 and 5 GHz, show a 
spectral flattening at v > 5 GHz. We identify 20 objects (19% of the sample) clearly 
characterized by an upturn spectrum, i.e., a spectrum falling at low frequencies (y < 
5 GHz) and inverted at higher frequencies. Spectra with high-frequency flattening or 
upturn shape are supposed to occur when the emission from the AGN compact core 
begins to dominate over the component from extended lobes. This picture fits well 
with the AGN unified scheme, for objects observed at intermediate viewing angles of 
the AGN jet. Finally, we discuss implications that this class of sources can have on 
future CMB observations at high resolution. 

Key words: 

Galaxies: active - quasars: general; Radio continuum: galaxies; Cosmic microwave 
background 



1 INTRODUCTION 

In recent years, extragalactic radio sources have received in- 
creased attention due to their contamination of the cosmic 
background radiation (CMB) temper ature fluctuations and 
polarization on s mall angular scales (|Toffolatti et al.ll2005l : 
iTucci et allbooih . At microwave wavelengths, this contri- 
bution is dominated by objects powered by active galactic 
nuclei (AGNs), such as quasars or BL Lacs. These sources 
are usually separated according to the slope of the spectrum 
at GHz frequencies (usually assumed a power law; hereafter 
S (x v a ): steep-spectrum if the spectral index a ^ —0.5 and 
flat-spectrum otherwise. The sources are thus divided into 
two classes, corresponding approximately to those where the 
radio flux is predominantly emitted from the extended lobes 
and t hose where th e compact core is the brightest compo- 
nent (|Peacockll 19851 ). 

Based on this scheme, evolutionary models for ra- 



dio sources provide accurate predictions for the num- 
ber counts at frequencies v < 8 GHz and down to flux 
densities of fractions of mjy ( Dunlop & Peacock 19901: 
i Toffolatti et al.lll998l: I Jackson fc Wal]||l999l : |Pe Zotti et all 
120051) . Recently, |Pe Zotti et all i|2005l ) updated and im- 
proved evolutionary models taking in to accounts new high - 
frequency data from the 15- GHz 9C (IWaldram et alj|200a ). 
18-GHz A TCA catalogues (TRicci et all 120041 ) and WMAP 
(first year; iBennett et alJl2003h . However, estimates of the 
radio source contribution to temperature fluctuations at 
cm/mm wavelengths, i.e. at wavelengths interesting for 
CMB observations, is still problematic. In fact, radio spec- 
tra of AGNs can be much more complex than a single power 
law, that may not provide a satisfactory description in a 
wide range of frequencies (see, e.g., simul taneous multi- 
freq uency measurement s of radio spectra by lTrushkinll2003l 
and iBolton et all |2004| ) . Different mechanisms can be re- 
sponsible for this: (i) a spectral steepening is expected due 
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to the more rapid energy loss of high-energy electrons with 
source age; (ii) a transition from optically thick to optically 
thin regimes can occur at high frequencies; (iii) at differ- 
ent wavelengths radio emission can be dominated by dif- 
ferent components characterised by a distinct spectral be- 
haviour. For example, AGNs with spectra peaked at GHz 
frequencies (GPS) are well known, while objects with spec- 
tra still inverted at frequencies as high as 20 -30 GHz are 
exp ected to fall at some peak frequency (see lO'Deal Il998l 
and |Pe Zotti et~aH 120051 for reviews on GPS and extreme 
GPS sources). 

Recently, indications of a class of radio sources with 
spectra that are "opposite " in respect to the GPS ones, are 
appearing in the literature. ISadler et al.l (|2006l ). using multi- 
frequency data of the ATCA 20-GHz pilot survey, observed 
objects with spectrum falling at low frequencies, but turn- 
ing up and begining to rise around 5 GHz, and called them 
"upturn" sources. T hey found that ab out 20% of the sources 
belongs to this class. iTrushkinl (|2003T ) identified in the first- 
year WMAP source catalogue 15 sources (corresponding to 
7%) with "combined spectra with a power low-frequency 
component and a flat (inverse) high-frequency component" . 
Upturn ob jects are observed also in recent obse rvations 
by ATCA fMassardi et aill2007l ; ISadler et al.ll2007h and by 
iPartridge et al.l ( 20071 ) in nearby galaxy clusters. 

The presence of a significant fraction of upturn sources 
in samples selected at few tens of GHz could have impor- 
tant implications for CMB observations and for models pre- 
dicting the radio source contribution at CMB frequencies. 
These sources, in fact, could appear in high-frequency sur- 
veys with flux densities much brighter than expected from 
low-frequency extrapolations. 

In this paper we carry out a multi-frequency spectral 
analysis of the sources detected by the exte nded version 
of th e Very Small Array (VSA) at 33 GHz (|Clearv et all 
2005, C05 hereafter). VSA provides a deep and nearly- 
complete survey of extragalactic radio sources with a flux 
limit of 20 mjy. A first discussion on the spectrum of these 
sources and their variability was carried out by C05, with 
an estimate of the source number counts at 33 GHz in the 
flux range 20-114 mjy. Here, we exploit surveys at differ- 
ent frequencies, as the NRAO VLA Sky Survey (NVSS, 
ICondon et all Il998l) at 1. 4 GHz, the Green Bank survey 



fGB6, iGregorv et al.lll996l) at 4.8 GHz and the Ryle Tele- 
scope CRT. IWaldram et all 120031 ) at 15 GHz, in order to 
study the spectrum of VSA radio sources in a quite extended 
range of frequencies. In particular, we focus our attention on 
sources whose spectrum shows a high-frequency flattening 
or an upturn behaviour, discussing a physical interpretation 
for such spectra. 



2 THE VSA SOURCE CATALOGUE 

The Very Small Array (VSA) consists of a 14-element inter- 
ferometer operating in the band 2 6-36 GHz, located at Teide 
Observatory in Tenerife (see, e.g. JWatson et ah| [2003). Two 
3.7 m antennas, spaced 9.2 m apart (called the "source sub- 
tractor"), are used to measure the flux density of sources 
present in the VSA fields. The sources to be followed by the 
subtractor are preselected from a catalogue at_15 GHz pro- 
vided by the Ryle Telescope (RT, see e.g. |jonesllT99ll and 
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Figure 1. Distribution of spectral indeces for steep— and flat- 
spectrum (defined between [1.4, 4.8] GHz) sources, calculated in 
frequency ranges [1.4, 4.8] GHz, [4.8, 15] GHz and [15, 33] GHz. 
Dotted histograms are excluding sources without 4.8-GHz mea- 
surements. Vertical dashed lines correspond to a = —0.5. 



IWaldram et al.l 120031 ) that surveyed the VSA fields in ad- 
vance of CM B observations. For the extended configurat ion 
of the VSA (|Grainge et alj 120031 ; iDickinson et al.1 \2004 ) it 
was required that all sources with flux S > 20 mjy were sub- 
tracted from CMB maps produced by the VSA main array. 
In order to satisfy this goal, the source subtractor followed 
all the sources detected by RT in the VSA fields up to the 
RT flux limit of 2 mjy. This ensures that all sources (except 
for extremely inverted ones) with S ^ 20 mjy at 33 GHz are 
observed by the source subtractor (C05). 

Our analysis is based on the catalogue of extragalactic 
radio sources presented by C05 and used to estimate the 
number counts at 33 GHz and to correct the source contri- 
bution to the VSA ex tended array power spectrum (see also 
IDickinson et al.ll2004l ) . Due to errors in the tex t of Table Al 
in C05 , we use the list of 33-GHz sources from lClearv et all 
(2008). The catalogue consists of 102 sources with flux den- 
sity between 20 and 114 mjy, within well-defined regions of 
the VSA fields that cover a total area of 0.044 sr. We con- 
sider also 3 further objects present in the VSA area and 
brigther than 114mJy. All these sources satisfy the pre- 
selection condition to be brighter than lOmJy at 15 GHz. 
As a consequence, the sample can not be considered strictly 
complete because there is the possibility that some sources 
with rising spectral index between 15 and 33 GHz (af| > 1) 
and S33 ^ 20 mjy exist in the VSA fields but it is not present 
in the sample. However, such inverted sources are rare and 
the loss of few of them would not affect the following anal- 
ysis. 
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Figure 2. "Two colour diagrams": the spectral index calculated at intermediate (left plot) and high (right plot) frequencies vs the value 
at low frequencies. The size of the circle is related to the flux densitiy at 15 GHz (left plot) and at 33 GHz (right plot). The vertical 
dashed line divides sources into steep— and flat-spectrum. Full dots in the right plot correspond to sources without 4.8 GHz flux; in this 
case, is substituted with o} 5 4 - 



3 SPECTRAL PROPERTIES IN THE 
FREQUENCY RANGE [1.4,33] GHZ 

Spectral properties of sources in the VSA sample can be 
studied using flux densities from surveys at lower frequen- 
cies: at 15 GHz from the RT; at 4.8 GHz from the GB6; and 
at 1.4 GHz from the NVSS. Only 11 objects do not have 
a counterpart in the GB6 survey, most of them character- 
ized by an inverted spectrum between 1.4 and 15 GHz and a 
low flux density at 1.4 GHz. We separate VSA sources into 
steep- and flat-spectrum sources according to their spec- 
tral index between 1.4 and 4.8 GHz (or 15 GHz in the case 
that 4.8-GHz flux density is not available). The sample is al- 
most equally divided into steep- and flat -spectrum sources, 
being 54 and 51 respectively, as predicted by evolutionary 
models for radio sou rces in this flux-density range (see, e.g., 
|Pe Zotti et~ai]|2005l ). 

Figure [l] shows the distribution of spectral indices for 
the two source populations in three different intervals of 
frequencies: [1.4, 4.8] GHz, [4.8, 15] GHz and [15, 33] GHz 
(hereafter indicated as low, intermediate and high frequen- 
cies respectively). We suppose = a\ 5 s = a} 5 4 for 
the 11 sources without measurements at 5 GHz. For flat- 
spectrum sources, the dispersion in the distribution of 
spectral indices increases at intermediate and high frequen- 
cies, with tails up to strongly inverted or steep indices 
(|q| > 1). However, the average value of a remains nearly 
constant with the frequency interval, 

(at!) = -0.05 (al%) = -0.15 (all) = -0.11 

Different remarks can be drawn for steep— spectrum 
sources. The distribution of a\'\ is very peaked between -1 
and -0.5, as expected because very steep-spectrum sources 
have low probability to be observed at high frequency. In the 



intermediate interval spectral indices spread both to positive 
and to very steep values, but keeping a similar average a. 
On the contrary, a clear "flattening" in spectra is observed 
at high frequencies, where more than 60% of sources have 
aff > —0.5 and 1/3 of steep sources has become inverted; 
the average value increases by +0.5 from the low- to the 
high-frequency interval: 

(atl) = -0.75 (al%) = -0.65 (all) = -0-25 

The presence of curvatures in source spectra is more 
evident in "two-colour diagrams": in Figure we compare 
the low-frequency spectral index, a\'%, with the one ob- 
tained at higher frequencies. In the left plot of Figure [2] 
we have a\ 5 s vs af;|: spectral indices seem to be distributed 
around the line a\% = af;|, with a bigger dispersion for 
flat-spectrum sources. Only in a few objects does the spec- 
tral index significantly change. No systematic trends can be 
observed and a single power law could statistically describe 
source spectra. This is not the case when we consider the 
high-frequency spectral index (aff vs af;|; see right plot). 
In most of the steep-spectrum sources the spectrum tends 
to flatten at high frequency or even to become inverted (i.e. 
with the typical upturn shape) . For flat-spectrum sources we 
notice a clear flattening only in objects with spectral index 
—0.5 < a\'X < — 0.3, that represent a transition from proper 
steep-spectrum (lobe-dominated AGN) and flat-spectrum 
(core-dominated AGN) sources. In general, if we consider 
sources with ct\% > —0.3, there is an average steepening in 
spectra at v > 5 GHz and very few objects keep the spec- 
trum flat or inverted at high frequencies. 

The radius of circles in Figure [2] is related to the flux 
density of sources: a visual inspection of the plots does not 
show any correlation between the position of sources in the 
two-colour plot and their flux density. To confirm this, we 
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Table 1. Number of sources with |Aa| = \au — cl\ > 0.3 (in 
brackets, > 0.5). 
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Figure 3. Distribution of spectral indices at low (left plot) and 
high (right plot) frequencies. Solid lines are the Gaussian fit of 
the distributions. 

calculate the linear correlation coefficient r between the flux 
density at 33 GHz, S33, and the spectral index variation 
Oil — cti'l- using all the sample, we find r = 0.027 with a 
propability for the null hypothesis of zero correlation of 0.76. 
This probability increases up to 0.8 (0.96) if we consider only 
sources with a\'X < (< —0.3). This is particularly impor- 
tant because the lack of correlation rules out that spectral 
flattening at high frequencies occurs only in faint objects, 
where uncertainties in flux densities are higher: as exam- 
ple, we find that the brightest source of the sample (source 
1734+3857) have a\l Z a l% > a tl- Moreover, a correla- 
tion between 33-GHz flux densities and spectral flattening 
should be expected if flattening is observed in objects that 
are just in the extreme of the spectral index distribution. 

3.1 Fit by two power law of source spectra 

For a better characterization of source spectra and their 
curvature (when present) , we suppose spectra to be well de- 
scribed by two different power laws. Let v m be the frequency 
at which the change in the spectral slope occurs. We define 
the low-frequency and the high-frequency (olh) spec- 

tral indices as the indices of the two best-fit power laws be- 
tween [1.4, Vm] GHz and [v m , 33] GHz respectively. The first 
step is to find out the frequency v m at which source spectra 
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Figure 4. Two colour diagram using spectral indices from two- 
power— law best fit. Open points are for objects whose ajj is calcu- 
lated between 15 and 33 GHz; solid points for objects whose ctjj 
is calculated between 4.8 and 33 GHz; crosses for objects with- 
out 4.8— GHz measurements. The diagonal dotted lines represent 
clh = &L + 0.3/0/— 0.3 respectively. 

are better fitted, considering that, in our case, v m can only 
be 4.8 or 15 GHz. Figure [3] shows the distribution of ai and 
olh of the VSA sources. The former presents a bimodal dis- 
tribution that can be fitted by two Gaussian distributions 
of mean —0.75 and —0.05 and dispersion 0.22 and 0.32 re- 
spectively. The two components can represent the two pop- 
ulations of AGN, steep- and flat-spectrum sources, clearly 
separated at low frequencies. On the contrary, at high fre- 
quencies the spectral index distribution is well fitted by only 
one Gaussian distribution with mean —0.25 and dispersion 
0.55, as result of the average flattening of steep-spectrum 
sources and the average steepening of flat-spectrum sources 
at high frequencies. 

The spectral properties of VSA sources between 1.4 and 
33 GHz can be now described by a single "two-colour dia- 
gram" using oll and olh (see Figure UJ. There we plot open 
(full) points for sources with v m = 15 (4.8) GHz; sources 
without 4.8-GHz measurements are indicated by a cross. 
Uncertainties in the spectral index are usually around 0.1- 
0.3 for both ol and an', however, they can increase more 
than 0.3 when olh is calculated only by 15 and 33-GHz data 
(open points), due to the uncertainties associated with the 
VSA measurements. 

The majority of sources (about 3/4; see Table [TJ are 
observed to change their spectral slope more than 0.3, indi- 
cating a possible curvature in the spectrum, while very few 
sources can be fitted by a single power law. An average flat- 
tening in steep-spectrum sources (ql ^ —0.5) is still con- 
firmed: 66% of these sources have Act = olh — oil ^ 0.3, and 
nearly 50% have Aa ^ 0.5. From Figure [4] we notice that 
the largest changes in a occur for sources with v m = 15 GHz 
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Figure 5. New GPS source candidates found in the VSA sample. 

(i.e. open points or crosses): these large differences could be 
overestimated due to the uncertainties associated with the 
VSA measurements and some caution is required. In the 
next section, we discuss about upturn sources and we deal 
with this issue. 

It is important to stress that few objects that are 
flat spectrum at GHz frequencies, keep their spectral slope 
flat or inverted up to 33 GHz in Figure [3] Instead, a gen- 
eral steepening is found. Looking at sources with cxl ~ 0, 
the typical steepening is about Aa ~ 0.3-0.5. This re- 
sult seems to indicate that, contrary to what normally 
assumed, a single power law with index a ~ is not, 
even statistically, a correct description of spectra for flat- 
spectrum sources at high frequencies. The contribution of 
flat-spectrum sources at the frequencies used for CMB obser- 
vations could thus be lo wer than predicted. Th is is in agree- 
ment with results from IWaldram et al. | (120071). with recent 
observations by th e ATCA telescope IjMassardi et al.l 120071; 
ISadler et al. 2007 ) and with spectra l properties of WMAP 
sources 1 Gonzalez-Nuevo et al.1120081) . 

Finally, we attempt to identify sources with spectra 
peaking at a frequency between 1.4 and 33 GHz, by the 
requirements oil > 0, ctn < and Aa < —0.3: there 
are 10 sources satisfying these conditions. Only 2 of them 
(0024+2911 and 1526+4201) are included in the 8 GPS 
sources identified by C05. Six of the remaining sources could 
be possible GPS candidates and their spectra are plotted in 
Figure [5] Further data are required to confirm their GPS 
shape. 



4 SOURCES WITH UPTURN SPECTRUM IN 
THE 33-GHZ VSA SAMPLE 

Having previously characterised the source spectra, in this 
section we look for objects with upturn spectra, or in other 
terms, with a minimum in the spectrum at v m = 4.8 or 



15 GHz. We say that a source has an upturn spec t rum if 
«l < and an ^ 0, in agreement with lSadler et al.l (2006), 
but with the additional condition an — oil > 0.5, in or- 
der to minimize spurious cases. Over a total of 105 sources, 
we find 27 objects with an upturn spectrum, correspond- 
ing to 2 6% of the sample (easily identified in Figure Q. 
While in ISadler et al.l (|2006h most of the objects classified 
as upturn are flat -spectrum at low frequencies with cxl ~ 0, 
in the VSA sample there are only 8 upturn sources with 
q_l > —0.5. This means that the majority of them are steep- 
spectrum sources at low frequencies that become flat or in- 
verted at v > 5 GHz, with a strong variation in the spectral 
index (in some cases > 1). 

The different survey resolutions ( < 1 arcmin for 
WENSS, NVSS and RT and ~3 arcmin for VSA and GB6) 
could affect the shape of spectra in the case of extended 
sources. Among the selected upturn sources, we find four 
that are resolved by NVSS, with the major axis between 44 
and 67 arcsec. Because the RT resolution is 25 arcsec, the 15- 
GHz flux density could be underestimated, even if integrated 
fluxes are considered. In particular, two extended sources, 
0010+3044 and 1240+5334, are characterized by a peculiar 
spectrum with a strong fall at 15 GHz (a 4 5 s < —1.5) that 
could be explained if power on angular scales higher than 
25 arcsec is missed by RT. Some caution is therefore needed 
for these objects and they are excluded from the analysis. 

It is important to rule out that upturn spectra could 
arise from bias in the 33-GHz flux density. It is known, in 
fact, that for flux-limited surveys the flux density of faint 
sources (i.e., close to the flux limit of the survey) could b e 
overestimated due to the Eddington bias (lEddingtonlll940h . 
Using Bayes' theorem, we can calculate the maximum likeli- 
hood estimator of the flux density for a source with detected 
flux d ensity So (|Hogg fc Turneriri998l : lLopez-Caniego et al.l 
l2007f ) as 



c So 

<->ML = — 



with r ^ 2a//3, 



(1) 



where r is signal-to-noise ratio of the detected flux density 
and (3 is the slope for the differential source count distribu- 
tion. For the VSA sample the estimated slope is f3 = 2.34 
(C05), implying that the Eq.[T]can be applied only to sources 
with signal-to-noise higher than 3 (this condition is veri- 
fied by ~ 97% of the VSA sources). In the few cases where 
S/N< 3, we consider 5ml = So/2. 

Now using the maximum-likelihood flux densities at 
33 GHz we calculate the new values for oll and an (the 
Eddington bias can be supposed negligible in the other sur- 
veys because of the high S/N ratio, at least for the VSA 
objects): after removing extended sources, 20 objects (19% 
of the sample) still verify the upturn conditions. We list 
them in Table [2] 

As discussed by C05, there is a deliberate deficiency of 
bright sources in the VSA fields. The first fields to be ob- 
served were selected to avoid sources predicted to be brighter 
than 250 mjy at 33 GHz; later fields used a more relaxed 
limit of 500 mjy. As a consequence, the fields are statisti- 
cally representative of the whole sky only between 20 and 46 
mjy for first fields and between 20 and 114 mjy for the later 
fields. Limiting to S ^ 46 mjy we find 17 upturn sources over 
a total of 86, corresponding to 20%, well in agreement with 
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Table 2. List of VSA sources with an upturn spectrum. The flux densities are at the frequencies: 325 MHz from the WENSS; 1.4 GHz 
from NVSS; 4.8 GHz from GB6; 15 GHz from RT; 33 GHz from VSA. The reported spectral indices are calculated using the VSA flux 
densities uncorrected by the selection bias. The last column reports the frequency of the minimum in spectra, as explained in the text. 



name S(mjy) ajj ajj — v, 





325 MHz 


1.4 GHz 


4.8 GHz 


15 GHz 


33 GHz 














GHz 


0013+2834 


— 


39.8±1.3 


23.8±3.7 


28.0±2. 


36.8±1.8 


-0. 


.42 


0. 


27 


0. 


.69 


4.8 


0023+2928 


85.±7. 


41.8±1.3 


26.2±5.2 


20.3±2. 


35.5±1.9 


-0. 


.31 


0. 


71 


1 


.02 


15 


0024+2724 


318. ±30. 


122.7±3.7 


41.4±4.1 


24.0±2. 


29.7±2.8 


-0. 


.73 


0. 


27 


1 


.00 


15 


0728+5325 


89.±7. 


60.±1.8 


36.8±4.1 


49.0±2. 


52.9±2.7 


-0. 


.10 


0. 


16 


0. 


.56 


4.8 


0728+5431 


72.±7. 


44.1±1.4 


40.0±4.6 


64.8±2. 


112.4±2.7 


-0. 


.08 


0. 


65 





.73 


4.8 


0931+3049 


111. ±8. 


67.6±2.1 


23.3±3.1 


24.0±2. 


32.6±6.7 


-0. 


.87 


0. 


10 


0. 


.97 


4.8 


0944+3115 


111. ±8. 


56.2±2.1 


39.5±5.1 


18.0±2. 


25.3±2.6 


-0. 


,16 


0. 


43 


0. 


.89 


15 


0946+3050 


62.±7. 


37.8±1.2 
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the percentage from the all sample. A smaller percentage is 
obtained in the area of 0.024 steradian, where there are 7 
upturn sources over 43 (16%) in the range 20-46 mjy and 2 
over 14 (14%) between 46 and 114mJy. 

In Figures [6HZl we show the spectra for the 20 upturn 
sources found in the VSA sample. For a visual inspection of 
spectra at low frequencies we include the flux densities at 
325 MHz as m easured by the W cstcrbork Northern Sky Sur- 
vey fWENSS. lRengelink et af1ll997D . We plot also the low- 
frequency best-fit power law: as expected for these sources, 
the extrapolation of the low-frequency power law to high 
frequencies generally underestimates the actual fluxes at 
33 GHz, in some cases up to a factor of 10, like for the 
objects 0931+3049 and 1229+5147. It is remarkable how, 
in general, the flux density at 325 MHz fits well the power 
law extrapolation from data at frequencies about one order 
of magnitude higher, confirming the reliability of the power 
law approximation for radio spectra at low frequencies. In 
few cases a flattening at 325 MHz is observed, due to a pos- 
sible effect of synchrotron self-absorpt ion or a low-energ y 
cutoff in the electron energy spectrum (jCarilli et al.lll99~D) . 

An interesting case is source 1734+3857, plotted in 
Figure [7] it is the brightest object in the VSA sample 
with a flux density a bit less than 1 Jy. This source has 
been detected in the WMAP 23, 33, 4 1 and 6 1-GHz chan- 
nels with a flux density ar ound 1.1 Jy |Hinshaw et al-lfeoOTl ; 
iLopez-Canieeo et al ] 120071 ). confirming the rising spectrum 
at v > 5 GHz observed by the VSA measurement. 

4.1 Variability in upturn sources 

Since observations of the surveys were carried out at differ- 
ent epochs, variability of radio sources could affect the shape 
of spectra. A question is therefore if the observed upturn 
spectra could be fictitiously generated by high variability in 



VSA sources. Two different cases can be discussed according 
to the value of u m . 

For upturn sources with spectral minimum at v m — 
15 GHz, only the VSA flux density does not agree with a 
spectral power-law fit. Because observations of the RT and 
VSA are separated by no more than two years, such sources 
should be high variable in a time scale of 1-2 years, and 
the VSA observations should correspond to a peak in the 
flux density. VSA observations cover a time scale from few 
months up to six months. As discussed by C05, the vari- 
ability of 72 sources could be examined: the majority of 
them (about 70%) varies by less than 25 per cent of the 
mean flux density, a 20% between 25 and 50 per cent and 
11% more than 50 per cent. Only three upturn sources with 
v m — 15 GHz are found in these 72 objects, and they show 
to have varied by less than 20 per cent. On the contrary, 
variability of at least 50 per cent is required to produce 
changes in the spectral index, calculated between 15 and 
33 GHz, bigger than 0.5, i.e. enough to explain upturn spec- 
tra. Moreover, high variable objects are usually associated 
to flat-spectrum sources (see the Figure 9 of C05 where a 
strong correlation between v ariability and spectr al index is 
found). This is confirmed bv lBolton et al.l (J2006) , who car- 
ried out a study of variability in a time scale between 1.5 
and 5 years: no steep-spectrum sources have been classified 
as variable whilst 50 per cent of the flat -spectrum sources 
have varied ( up to a 70% of the mean flux density). Also 
ISadler et all (|2006f ) found that at 20 GHz the general level 
of variability is quite low in a one-year interval, with only 
five sources varied by more than 30%. To conclude, it seems 
quite unlikely that variability could explain upturn spectra 
with minimum at 15 GHz. 

On the other hand, the total set of data covers nearly 
twenty years (the GB6 survey was taken in 1986-87, WENSS 
in 1991-93 and NVSS between 1993 and 1996, RT in 2000- 
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Figure 6. Sources with upturn spectra. Dotted lines are the power laws that best fit spectra at low frequencies. 



2001 while VSA in 2002-03). Variability on time scales of 
years could be relevant for sources with spectra upturning at 
Vm = 4.8 GHz: because there are no observations following 
these sources on time scales of ten years or more, we are not 
able to exclude a significant long-scale variability in their 
flux density. We consider the "worst" case, in which the flux 
densities measured by GB6 at 4.8 GHz (the most distant in 
time from VSA measurements) correspond to a minimum 
in the source emission, giving a steeper low-frequency spec- 
tral index. Ignoring the 4.8-GHz fluxes, we take ol = u{ 5 4 : 
over the 10 upturn sources with the minimum at 4.8 GHz 
only 2 still satisfied our upturning conditions, even if 7 of 
them keep an upturn shape with the less strict conditions 
ql < 0, an > and Aq > 0.3. Therefore, in principle, a 
long-scale variability could explain the spectral curvature 



observed in sources with minimum at v m = 4.8 GHz (but 
not in ones with v m = 15 GHz) requiring a large fractional 
variability, in many cases of the order of 100 per cent at 
4.8 GHz. Because at the moment there are no observational 
data supporting this hypothesis and because variability is 
not able to explain all the upturn spectra found in the VSA 
sample, we are still confident that the observed upturn spec- 
tra represent intrinsic properties of radio sources. Finally, 
we observe that in most of upturn sources the WENSS flux 
densities at 325 MHz are in very well agreement with extrap- 
olations based on the estimated a l , indicating no variability 
on a time scale of years at these frequencies. 
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Figu re 7. Sources with upturn spectra (as in Figure [B}. For source 1734+3857 we plot also the flux densities obtained by 
iLopez- Canicgo et al. (2007) using the 3— year WMAP data (open points). 



5 DISCUSSION ON FLATTENING/UPTURN 
SOURCES 

Although not explicitly considered, sources with flattening 
or upturn spectra are co nsistent with the unif i ed scheme for 
rad io-loud AGNs (e.g., lOrr fc Browne! Il982l . iBarthell Il989l 
an d lUrrv fc Padovanilll995l )T According to this picture, spec- 
tral properties of the AGN radio emission depend on the 
degree of the alignment of relativistic jets with the line of 
sight (hereafter the "viewing angle"): in a simplistic way, 
steep spectra are observed in misaligned objects, where ra- 
dio emission is dominated by extended optically-thin lobes; 
on the other hand, if the line of sight is enough close to the 
AGN jet, the optically-thick and beamed emission of the 
AGN compact core emerges with the typical flat spectrum 



|Peacock|[l985h . The critical viewing angle 6 C that defines 
the division between flat- and steep-spectrum objects at 
GHz frequencies is us ually considered around 10° for quasars 
and 3 0° for BL Lacs (|Ghisellini et al.lll993l ; Ijackson fc Walll 
Il999l ). Therefore, for angles 8 not much larger than 9 C , we ex- 
pect a transition from steep to flat spectra at some frequency 
v > 5 GHz, when the compact emission begins to dominate 
over the exten ded one (whose spectr um falls rapidly with 
the frequency). I Jackson fc Wa ll (1999) employed this idea in 
order to extrapolate the number counts of radio sources at 
151 MHz to 5 GHz, introducing the "beamed" contribution 
of parent populations of blazars. This picture is in agree- 
ment with high-resolution observations: AGN radio struc- 
tures are commonly resolved into three basic components 
(nuclear core, large scale jets and outer lobes), while ex- 
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tended r adio structures without a c ontribution from a core 
are rare (IPearson fc Rcadhcad 1988). Moreover, the relative 
importance of the AGN core emissions is obse rved to in- 
creas e at higher and high er frequencies (see, e.g.. lBach et al.l 
|2004| for Cygnus A and iKharb et al.l I2OO7I for a sample of 
Fanaroff- Riley II) . 

Now we discuss more quantitatively how the unified 
scheme applies to our case. We consider the observed flux 
density of VSA sources to be the sum of two compo- 
nents, an extended steep-spectrum and a flat -spectrum 
core-dominated emission. For simplicity, we suppose that 
the VSA sample is composed of quasars, being the domi- 
nant population at that frequency and flux range. Let R° bs 
be the observed core-to-extended flux ratio at 5 GHz. At 
a frequency v > 5 GHz the core component dominates the 
total emission if 



R 



obs 



V Sext I 



> (sis) < 2 » 

where a ext and a c are the spectral index of the two com- 
ponents, extended and compact. In the following, we adopt 
as spectral indices the typical values for steep- and flat- 
spectrum s ources at low freque ncies, a ex t = —0.75 and 
q c = -0.1 (|De Zotti et al.ll2005h . The value for a ext could 
be conservative at frequencies v > 10 GHz because it does 
not take into account possible steepening due to the electron 
ageing. 

Objects with flattening or upturn spectrum in the 
VSA sample should be characterized by S c /S ex t > 1 at 
33 GHz, corresponding to R° bs > 0.3. A relevant fraction 
of them, moreover, shows flattening/inverted spectra al- 
ready at 15 GHz, requiring R° bs > 0.5. If the core emis- 
sion comes out from a relativistic jet of bulk velocity /3c, 
the observed flux density is Doppler-enhanced by 8 v ~ a " , 
where 8 — [7(1 — (3cos6)]~ is the Doppler factor and 
7 = (1 — /3 2 ) -1 / 2 is the Lorentz factor. The parameter 
p depends on the emission model and is equal to 3 for a 
moving, isotropic sou rce and 2 for a continuou s jet (see 
iGhisellini et~ai1 Il993l and lUrrv fc Padovanil Il995l for a de- 
tailed discussion). Then, the relation between the observed 
(R° bs ) and the intrinsic-unbeamed (/) core-to-lobe flux ra- 
tio for a source at redshift z is given, after the K-correction, 
by 



R obs = / 8 p ~ a < (1 + . 



(•3) 



In the unified scheme, the value of / is the same for quasars 
and for their parent population, the FRII galaxies. Because 
there is no evidence that the parameter / correlates with the 
radio power of FRII galaxies l| Jackson fc WalCl 999). we can 
use t he best-fit FRII value f = 5 x 10~ 3 (|Urrv fc Padovanil 
1 19951 ) for both steep- and flat-spectrum quasars, indepen- 
dently of their luminosity. 

Using Eq. [3] now we are able to provide an indicative 
estimate of R° bs for quasars observed at a viewing angle 
10° ^ 9 20°: supposing a Lorentz factor between 2 and 
15, we obtain 
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The K-correction is of little relevance, being a factor ~ 1.6 
for 2 = 1 (the redshift we use). Fixing the viewing angle and 
for the considered parameters, the value of R° bs scatters 
over a range of a factor 30. In any case, we note that for 
9 = 10° R° bs is compatible with the values required to have 
an upturn spectrum, independently of the 7 and p used. 
This is not the case when 9 is near to 20° (which should be 
considered as upper limit) and upturn spectra are possible 
only for low values of the Lorentz factor (7 < 5) and p = 3. 
The case of 9 = 15° may be more conservative because the 
condition that R° bs = 0.3 is verified independently of p, at 
least for 7 < 10. Therefore, according to this model and 
for the range of frequencies investigated, a viewing angle 
of 9 < 15° is required to observe quasars with a flattening 
spectrum. 



5.1 Predicted number of flattening/upturn 
sources in the VSA sample 

A prediction of how many upturn sources should be ex- 
pected in a sample selected at 33 GHz is quite difficult, not 
only for the scattering in the physical parameters used in 
Eq. [3] but also for selection effects in flux-limited samples. 
However, following the previous discussion, a qualitative es- 
timate is possible. We want to stress that nearly all the 
VSA upturn sources (18/20) have ai < —0.3. Sources with 
—0.5 < o?l < —0.3, although not strictly steep spectrum ac- 
cording the usual definition, can be considered a transition 
between lobe- and core-dominated quasars. In the following 
we do not make a distinction between them. 

We suppose that all the sources with steep spectra at 
low frequencies observed at viewing angles 10° ^ 9 ^ 15° 
show a flattening or rising spectrum at v > 5 GHz. If lobe- 
dominated quasars are in general observed at 10° ^ 9 ^ 40° 
l|Ghisellini et al.lll993l ). we should find that about 1/6 of 
steep-spectrum sources have a flattening or upturn spec- 
trum at high frequencies (in the hypothesis that the proba- 
bility to observe a radio quasars is independent of 9). This 
estimate should be only applied to samples selected at few 
GHz. In flux-limited samples selected at high frequencies, 
this percentage could be higher because of missing very steep 
spectrum sources. As example, in the VSA sample we have 
found 18 upturn objects out of a total of 69 sources which 
are steep {at < —0.3) at low frequencies. 

In order to establish if the number of flattening/upturn 
sources expected from the model is in agrement with the 
number observed in the VSA sample, we proceed in the fol- 
lowing way: 

• we suppose that for flattening/upturn sources the aver- 
age "high-frequency" spectral index is (off) = 0. 

• Using the NVSS and GB6 surveys we select inside the 
VSA fields all sources with a\'\ < —0.3 and flux density 
higher than 20mJy at 4.8 GHz. We find 370 objects. 

• Extrapolating the flux density to 15 GHz using a\'%, 
we find that 149 of these 370 sources have S15 20mJy. 
According to the above discussion, we expect that 1/6 of 
them (i.e. 25 objects) have flat or inverted spectrum at high 
frequencies {an ~ 0), and are observed in the VSA sample. 

• Moreover, the rest of the 221 steep-spectrum sources 
(with S4.8 ^ 20mJy but the extrapolated $15 < 20mJy) 
can also be observed in the VSA sample if their spectra are 
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flat starting from v — 4.8 GHz. These objects correspond to 
the previous cases with v m = 4.8 GHz, and they are more 
or less one half of the total number of flattening/upturn 
sources. Therefore, considering a fraction of 1/12, we have 
a further contribution of 18 objects. 

To resume, the number of sources with flatten- 
ing/upturn spectrum predicted by our estimates in the VSA 
area at 33 GHz is 43. This value should be compared with 
the number of steep-spectrum (at < —0.3) sources with 
Aa = an — ctL ^ 0.3 found in the VSA sample. This dif- 
ference is in fact approximately the change of the spectral 
index in the case S c > S ex t at 33 GHz (using the a c and 
a ext assumed above). As seen in Table[T] there are 43 sources 
with these characteristics, as from our predictions. A smaller 
number, 37, is found if we correct the 33-GHz VSA flux den- 
sities by the Eddington bias. 

If we want to limit our predictions only to sources 
with an upturn behaviour, i.e. rising at high frequencies, 
we should require that a c > and S c S> S ex t at 33 GHz. 
The fraction of flat-spectrum sources with a > can be ob- 
tained from samples of sour ces whose emission is dominated 
by the AGN compact core: iRicci et al.l l|2006l ) showed that 
for a complete sample of flat -spectrum sources at 5 GHz, the 
distribution of the spectral index calculated between 2.7 and 
5 GHz is well described by a Gaussian with (a) = —0.1 and 
a a = 0.3. These values are similar to parameters describing 
the spectral index distribution at low frequencies for flat- 
spectrum sources in the VSA sample (see Figure |3J. Using 
these parameters in a Gaussian distribution, we find that 
about 38 per cent of spectral indices are higher than 0. This 
means that 16 of the 43 predicted objects with flattening 
spectrum is expected to be upturn between 5 and 33 GHz, 
again in agreement with the results from the VSA sample. 



6 UPTURN SOURCES FROM 

HIGH-FREQUENCY SURVEYS 

In this section we search for sources with upturn spectra 
in other surveys at frequencies comparable or higher than 
the VSA one. The fraction of "upturn" sources is expected 
to depend on th e selection frequency an d the flux-density 
range of surveys l|Kellermann et al"1ll968T ). A direct compar- 
ison with the VSA results should take into account this. 

• WMAP source catalogue. The WMAP mission has 
produced the first all-sky surveys of extragalactic sources 
at 23, 33, 41, 61 and 94 GHz and at Jy flux densities 
dBennett et al 1 l2003l ; iHinshaw et al]|2007l ). In this interval 
of frequencies, spectral properties of extragalactic sources 
are still poorly inves tigated. From 3-years WMAP maps, 
IHinshaw et al.l l|2007i ) provided a catalogue of 323 objects 
ab ove a flux limit of abou t 1 Jy. This catalogue was improved 
by lLopez-Caniego et al.l l|2007i ) using a non-blind method 
based on the second member of the Mexican Hat Wavelet 
Family. They found 368 sources with a 5-er detection in at 
least one WMAP channel and with a completeness flux limit 
of around 1.1 Jy. 

We complement the WMAP flux densities by measure- 
ments from low-frequency surveys: in the northern hemi- 
sphere from the NVSS and GB6 surveys; in th e southern 
hemisphere from the SUMSS (|Mauch at al.ll2003T l and PMN 



l| Wright et al.ll 19941 : IWright et ai1ll996T l surveys at 843 MHz 
and 4.85 GHz respectively. We use these data to define the 
low-frequency spectral index (q.l) for WMAP sources. Ob- 
jects without a counterpart both in NVSS/SUMSS and in 
GB6/PMN or without a clear identification have been dis- 
carded from the analysis: 248 objects remain, 114 with 
low-frequency spectral index ql < and only 22 with 
q_l —0.5. If the high-frequency spectral index is calcu- 
lated between 4.8 GHz and the lowest frequency at which 
a source has been detected by WMAP, we find 25 upturn 
sources (5 with «l ^ —0.5) according to the previous def- 
inition. Some examples of th ese with counterpa rts in the 
multi-frequency cat alogue bv | Kuehr et al.l lll98l[) or in the 
ATCA observations l|Ricci et alj|2006l : iMassardi et al.ll2007f ) 
are shown in Figure [5] 

The detection of objects with upturn spectrum in the 
WMAP catalogue is an important confirmation of the pres- 
ence of this class of sources also at Jy flux levels. Their 
fraction is about 10%, lower than in the VSA sample. This 
is not unexpected due to the different range of flux densi- 
ties considered by the two experiments. At Jy level, in fact, 
the source population is dominated by flat-spectrum sources 
jDe Zotti et al l I2005D and the fraction of objects with in- 
verted spectrum at low frequency is particularly high. In the 
WMAP catalogue more than 50% of sources have ol > 0, 
compared to 17% in the VSA sample. 

For a better comparison between the two samples we con- 
sider only the WMAP sources with a 5-a detection in the 
33-GHz channel, in order to achieve a sample selected at 
the same frequency of the VSA. Then, we compute ql and 
an as in section 3.1 (using the 23-GHz WMAP channel 
instead of the 15-GHz frequency) . Now, if we limit the com- 
parison only to objects with decreasing spectrum at low fre- 
quency (ol < 0), we find about 20% of upturn sources in 
the WMAP catalogue, very close to 23% in the VSA sample. 

• The 20 GHz ATCA Survey. A pilot survey of 
the Australia Tele scope 20 GHz survey was presented by 
ISadler et al.l l|2006h . consisting of a flux-limited sample (S ^ 
100 mjy) of 173 radio sources selected at 20 GHz, includ- 
ing near-simultaneous flux measurements at 4.8, 8.6 and 
18 GHz. In the paper, usin g also the flux den sities provided 
by the SUMSS catalogue l|Mauch at al.ll2003h , they carried 
out an accurate analysis of source spectra in the frequency 
range 0.8-20 GHz: about 20% of source spectra were clas- 
sified upturn (defined as Qq.s < and ag 8 > 0). However, 
introducing our stricter condition that ag 8 — ao.s ^ 0-5, the 
fraction is reduced to 14%. 

Recently, a complete sample of bright sources (S20GHZ > 
0.5 Jy), carrie d out with ATCA, was presented by 
IMassardi etaH (|2007l ). For 218 of them, near simultane- 
ous observations at 4 .8, 8.6 and 20 GHz were also provided. 
IMassardi et al.l (|2007T ) studied the source spectral behaviour 
in the frequency range 5-20 GHz: a general steepening in 
spectra is observed, while very few upturn sources are found 
(only three according our definition; see their Figure 4). Ex- 
tending the analysis to the NVSS and SUMSS measurements 
at 1.4 and 0.8 GHz, the number of sources with upturn spec- 
tra increases to 19 (8% of the sample), indicating that their 
spectral minimum is at < 5 GHz. However, as for the WMAP 
case, if we exclude objects with inverted spectrum at low fre- 
quency (that are about 60% of the sample) , the fraction of 
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Figure 8. Example s of upturn spectra in the WMAP catalogue. 
Counterparts in the iKuehr et al,l ill 98 if ) and in the ATCA sam- 
ples are shown as cross points and open squares respectively. Full 
points are from WMAP and from SUMSS/NVSS and PMN/GB6 
at low frequencies. 



upturn sources is of the order of 20%, well in agreement with 
the VSA and WMAP one. 

• R adio follow— up of the 9C survey. iBolton et al.l 
l|2004h carried out a follow-up of 176 sources from the 
15-GHz 9th Cambridge survey in the frequency range 
[1.4, 43] GHz. In particular, they made simultaneous mea- 
surements at 1.4, 4.8, 22 and 43 GHz with the VLA and at 
15 GHz with t he Ryle Telescope. A s discussed by the au- 
thors (see also lWaldram et al.l 120071 ). there is a clear trend 
for spectra to be steeper at the high frequencies: an average 
steepening is already observed at 22 GHz, but it becomes 
especially strong in the range 15-43 GHz (the median spec- 
tral index is -0.89), with a large number of objects with 
extremely steep spectra, less than —1.5. We find only 22 
(14%) inverted sources between 4.8 and 22 GHz and very 
few objects inverted up to 43 GHz or upturn. 

This strong steepening at v > 15 GHz has not been ob- 
served in other surveys (ATCA measurements found a steep- 
ening in 20-95 GHz spectral indices, however the median 
values is -0.39: see ISadler et al] l2007h and it is clearly in 
disagreement with VSA results. Such a discrepancy may be 
related to the different resolution of the instruments, e.g. a 
few arcsec for the VLA at 22-43 GHz, 25 arcsec for the Ryle 
Telescope and just 3 arcmin for the VSA source subtractor. 



7 DISCUSSION AND CONCLUSIONS 

An accurate control of the contribution of extragalactic ra- 
dio sources at microwave wavelengths will be extremely im- 
portant for future CMB experiments at arcminute resolution 
like Planck and for experiments aim ed to the detection of the 
cosmological B-mode polarization (|Tucci et a 1. 1 [201)5;. The 



VSA source-subtractor observations provide a very useful 
nearly-complete catalogue of "faint" (20-100 mjy) extra- 
galactic radio sources at a frequency used for CMB mea- 
surements. In this paper, using additional measurements at 
lower frequencies, we investigate the spectral properties of 
the VSA catalogue in the frequency range 1.4-33 GHz. 

The most interesting result we find is related to 
the high-frequency spectral behaviour of sources whose 
spectrum is steep at GHz frequencies. Typical models 
for this source pop ulation (see, e.g., the recent work by 
|Pe ZottTe t al. 2005]) predict a spectral steepening at high 
frequencies due to the rapid ageing of high-energy electrons, 
responsible of the emission at v 3> 1 GHz. On the contrary, 
from the 33-GHz VSA catalogue we observe that the spec- 
trum of the majority of them flattens or even becomes in- 
verted at v > 5 GHz. In particular, in about 50 per cent of 
steep-spectrum sources the spectral flattening at high fre- 
quencies is more than 0.5. Moreover, we find objects with 
upturn spectra, i.e. steep at low frequencies but rising be- 
tween 5 < v < 33 GHz: 20 sources clearly show such spectral 
behaviour, corresponding to 19 per cent of the total sample. 

The large fraction of sources characterized by spectral 
flattening and the lack of correlation between flux density 
and spectral curvature are indications that we are observ- 
ing intrinsic properties of radio sources, and not the result 
of selection effects or bias. Variability can affect the shape 
of spectra observed in measurements that are not simulta- 
neous. However, it is quite unlikely that variability is re- 
sponsible of the average spectral flattening found in VSA 
steep-spectrum sources. In fact, it would require large vari- 
ations in flux densities on time scales of both 1-2 years and a 
decade, for a source popul ation usually observed to exhibit 
little variability (see, e.g., I Sadler et al. I l2006l . IBolton et al.l 

l20oet) . 

In the paper, we show that source spectra with flatten- 
ing or upturn behaviour are completely consistent with the 
"unified model" for AGNs: steep spectra in radio sources 
can be physically understood as emission from extended 
optically-thin radio lobes, while a flat -spectrum component 
arises from the central compact and optically-thick core 
of radio sources. Due to the rapid decrease in intensity of 
the diffuse emission with frequency, the emission from the 
compact core of AGNs may start to dominate at frequency 
v 3> 1 GHz if the radio jet ax is lies enough close to the line of 
sight (Ja ckson fc Walllll999l ). According to the typical phys- 
ical parameters for quasars, it can occur for v > 5 GHz when 
the viewing angle is 6 < 15°. 

Finally, VSA data confirm that a simple power-law ex- 
trapolation of low-frequency flux densities to mm wave- 
lengths is not reliable: this is not only because a fraction 
of sources predicted at high frequencies are not found, but 
also because a significant number of sources that would not 
be expected are detected. For example, using a simple ex- 
trapolation of the spectral slope between 1.4 and 4.8 GHz 
to 33 GHz, we find that 34 sources present in the VSA cat- 
alogue would be predicted to be fainter than the flux limit 
of 20mJy and not included in the sample. From surveys in 
literature, the percentage of these "unexpected" objects is 
observed to steadily increas e with the frequenc y. While it is 



roughly 10 % at 15 GHz (seelTavlor et al.ll200lh and 18% at 



20 GHz (see lSadler et al.1 12006), this percentage increases up 
to 32% in the VSA sample. 
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We conclude that radio sources with flattening or up- 
turn spectrum could be an important contaminant to high- 
resolution CMB observations providing to an extra contribu- 
tion to source number counts respect to predictions based 
on surveys at 1-5 GHz. At 33 GHz this extra contribution 
is not negligible: limited to the 20 upturn sources found in 
the VSA sample, we substitute their measured flux density 
with the one obtained from a low-frequency extrapolation; 
in this case, the VSA number counts at 33 GHz would keep 
the same slope as found by C05 but rescaled down by 15 per 
cent. This is for flux densities 20 ^ S < 100 mjy. At higher 
fluxes the fraction of upturn sources is observed to be lower, 
around 10 per cent in the WMAP and ATCA source cata- 
logues at Jy levels. If confirmed, it would suggest an extra 
contribute from upturn sources of about 10-15% to temper- 
ature fluctuations generated by extragalactic radio sources 
at 33 GHz. At higher frequencies the relevance of this class 
of sources is expected to increase. 

In the near future a number of new ground-based sur- 
veys at 20-30 GHz will be available, as, for example, the 
super-extended version of VSA, the Australia Telescope 
20 GHz survey and the One Centimetre Receiver Array 
(OCRA) project. They will be able to more deeply inves- 
tigate radio source spectra in large intervals of fluxes and to 
confirm the results discussed in this paper. However, obser- 
vations at frequencies higher than 30 GHz are required for a 
better characterization of upturn sources. In this sense, the 
ESA's Planck satellite will be extremely important, carry- 
ing out all-sky surveys from 30 to 860 GHz with a sensitiv- 
ity several times better than WMAP, and extending high- 
frequency radio surveys to fluxes of few hundred of mjy. 
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